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ABSTRACT 

The Gould Belt Legacy Survey will survey nearby star- forming regions (within 500 pc), us- 
ing HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common- 
User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope 
(JCMT). This paper describes the initial data obtained using HARP to observe 12 CO, 13 CO 
and C 18 / = 3 -> 2 towards two regions in Orion B, NGC 2024 and NGC 2071. We de- 
scribe the physical characteristics of the two clouds, calculating temperatures and opacities 
utilizing all three isotopologues. We find good agreement between temperatures calculated 
from CO and from dust emission in the dense, energetic regions. We determine the mass 
and energetics of the clouds, and of the high-velocity material seen in 12 CO emission, and 
compare the relative energetics of the high- and low- velocity material in the two clouds. We 
present a CLUMPFIND analysis of the 13 CO condensations. The slope of the condensation 
mass functions, at the high-mass ends, is similar to the slope of the initial mass function. 

Key words: stars: formation - ISM:kinematics and dynamics - submillimetre - molecular 
data 

1 INTRODUCTION 
1.1 The Gould Belt Legacy Survey 

The Gould Belt Legacy Survey (GBS, | Ward-Thompson et al.| 
2007 ) has been awarded 612 hours of time on the James Clerk 
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Maxwell Telescope (JCMT) to survey nearby star-forming regions 
(within 500 pc), using HARP (Heterodyne Array Receiver Pro- 
gramme, [Buckle£tal]200^ Common- 
User Bolo meter Array 2, |Holland et al.||2006] ) and POL-2 (Po- 
larimeter 2, Bastien, Jenness & Molnar 2005). The HARP com- 
ponent will observe a large typical sample of prestellar and pro- 
tostellar sources in three CO isotopologues. CO J = 3 and the 
same transition in 13 CO and C 18 are excited in the physical condi- 
tions of star-forming molecular clouds, where typical temperatures 
and densities are 10-50 K and 10 4 -10 5 cm 5 . iZ CO observations 
are used to search for and map any high- velocity outflows present, 
while the isotopologues measure the line widths and velocity pro- 
files in the cores and filaments, resolving detailed kinematic and 
density properties of the cores. The GBS aims to provide a large 
and unbiased sample of star-forming material in the solar vicinity at 
relatively high spatial (8-14 arcsec) and spectral (0.05-1.0 kms -1 ) 
resolution. The GBS will trace the very earliest stages of star forma- 
tion, through submillimetre continuum imaging using SCUBA-2 
observations of low-temperature, high-column density regions, and 
provide an inventory of all the protostellar objects contained in the 
nearby molecular clouds of the Gould Belt, covering ~ 700 deg 2 . 
The HARP observations trace the kinematics of molecular gas in 
the cores and clusters at the same spatial resolution as the 850 \im 
SCUBA-2 observations of the dust emission. The key science goals 
of the HARP component of the GBS are ( Ward-Thompson Ik al.| 
|2007] >: 

• To search for and map high- velocity outflows in the cores in 
order to differentiate between starless and protostellar cores. 

• To derive simple constraints on the column density and CO 
depletion in these cores. 

• To investigate support mechanisms and cloud/cluster/core 
evolution. 

• To characterize the cloud kinematics in a wide variety of en- 
vironments and investigate the role and evolution of turbulence in 
star formation. 

The HARP targets for the GBS are a sample of the cloud re- 
gions that will be observed by SCUBA-2, with sizes ~0.8 deg 2 , 
and containing filaments and clustered star formation. This paper 
describes the initial HARP data obtained for two targets in Orion 
B, NGC 2024 and NGC 2071 . 



1.2 Orion B 

The Orion B cloud complex is the closest region of high-mass star 
formation, lying at a distance of 415 pc ( Anthony-Twarog 1982, 
Ment en et al.|20 07 ). Five main star formation regions were identi- 
fied from early wide field CO surveys: NCG 2024, NGC 2071, LBS 
23 (HH 24), NGC 2068 and NGC 2023 ([Tucker, Kutner & Thad- 
deus|1973l[Kutner et al.|l 977)1 White & Phillips|1981||Wilson et al. 
2005j r The cloud complex subtends an area of ~ 26 deg 2 , which 
corresponds to ~ 1.5 kpc 2 , and has a mass of - 8 x 10 4 M . 
The main star forming regions are clearly visible in the unbiased 



surveys for young stella r objects (|Lada, Bally & S tark 1991 [|Lis,| 
[ Menten & Zylkal[T999l fM otte et al.||2001[ |Mitchell et al.|(200Tj 
|johnstone, Matthews & Mitchell||2006| hereafter J06) and in the 
dense gas tracer CS ( |Lada et al.|1991| >. Observations of individual 
sources in CO lines revealed a complex structure of outflows and 
gas m otions ( [White, Phillips & Watt|19~8T||Phillips, White & Ade] 
1982), which are now known to be associated with extensive star 
formation activity (Flaherty & Muzerolle 2008). A detailed sum- 
mary of the Orion B regions can be found in the Handbook of Star 
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Figure 1. NGC 2024 and NGC 2071 shown within the Orion B cloud com- 
plex seen in visual extinction, A v from the atlas of Dobashi et al. ( 2005). 
The contours are drawn at A v =l,2,3,4,5,6 mag. 



Forming Regions (see chapters by |Gibb|2009| |Bally|2009l [Meyer 
|2009t . Th e two GBS Orion B targets, NGC 2024 and NGC 2071, 
for which we have obtained the first data are shown as boxes in 
Fig. ^overlaid on an extinction map of the Orion B cloud complex. 



1.2.1 NGC 2024 



NGC 2024 (The Flame Nebula, Anthon y-Twarog|1 982 ) is a bright 
emission nebula that is crossed by a prominent dust lane. Intense 
radio continuum and recombination line emission ( jKriigel et al, 
T982) |Barnes et al.|[T989| |Bik et al.|[2003| |Rodriguez, Gomez & 



Reipurth 2003), plus the detection of an infrared cluster, suggested 
the presence of ionizing star(s) of spectral type in the range 08 V- 
B2 V ([LadaT Bally & Stark|1991[|Comer6n, Rieke, & Rieke 



1996 



2003 



IGiannini et aTTT oOO ; Hais ch, Lada & Lada|20001 |Bik et al. 
|Kandori et al.|2007| ). 

Millimetre and submillimetre wavelength observations re- 
vealed a number of compact condensations that were believed to 
mark star forming sites. NGC 2024 has also been shown to con- 
tain a number of protostars positioned along a ridge of enhanced 
star formation activity, and its associated dense molecular ma- 
terial distributed as a filamentary structure ( [Mezger et al.| 1988 
1992, Visser et al. 1998 ). Far-infrared and submillimetre polari- 
metric observations tracing the magnetic field structure of the star- 
forming ridge suggest it contains helical magnetic fields surround- 
ing a curved filamentary cloud, possibly driven by the ionization 



© 2009 RAS, MNRAS 000,[l]j20] 



GBS HARP Survey: Orion B 3 



front of the expanding Hll region (H ildebrand et a l.|1995[[D"o tson 
et al.|2000[|Greaves, Holland & Ward-Thompson|2001[|Matthews,| 
Fiege, & Moriarty-Schieven 2002). Previous molecular line obser- 
vations of NGC 2024 have included studies of the outflow and ther- 
mal properties of th e gas (|Watt et al.| 1979 1| White & Phillips|1981| 
Philli ps et al.|1988||Bames et al.|19891|Richer et al.|1989[|Barnes| 
& CrutcherH 19901 IRicher, Hills & Padm an|| 1 992) |C hernin|| 19961 
Mangum, Wootten & Barsony 1999, Snell et al. 2000), revealing a 
complex and boisterous medium permeated with outflows, turbu- 
lence and complex gas dynamics. 

1.2.2 NGC 2071 

NGC 2071 is associated with a reflection nebula, molecular cloud 
and infrared cluster NGC 207 1IR (White & Phill ips| 1981] [W hite 
let al.|1981l|Persson et al.|1981[|Motte et al.|2001H Johnstone etall 
2001 hereafter J01). The main infrared cluster has been resolved 
into eight distinct near-infrared sources having a total luminosity 



of 520 L ( [Burner et al.|[T990l |Walther et al.p993) . The bipolar 
molecular outflow from NGC 207 1 has been extensively studied in 
CO emission ([White & PhTTTipsl[T98T] |Snell et al.p9Pl |Scoville| 



et al. 


|1986[ |Moriarty-Schieven, Hughes & Snell||1989| |Kitamura 


et al. 


1990| |Houde et al.|2001 , Eisloffel 2000), and is found to be 


amongst the most energetic bipolar outflows currently known, and 
extends for ~ 2 pc (Margulis & Snell 1989, Chernin & Masson 


1992 


|Chernin & Welch 1995, Stojimirovic, Snell & Narayanan 


2008 


i. This outflow is also seen in the shock excited 2.12 |im H2 


line ( 


Persson et al.|1981||Lane & Bally|1986||Burton, Geballe & 


Brand 1989, Garden, Russell & Burton 1990, Aspin, Sandell & 



Walther 1992). 

Snell & Bally (1986) have reported the presence of several ex- 
tended radio continuum sources. Subsequent observations of H2O 



masers and the associated radio continuum emission ( [Smith &[ 
|Beck|1994[|Torrelles et al.|l998i |Seth, Greenhill & Holder|2002| > 
suggest that the infrared sources IRS 1 and IRS 3 have strong 
winds, and are surrounded by rotating circumstellar disks. Fur- 
ther outflow activity around IRS 3 was suggested by its H2 emis- 
sion ( [Aspin et al.|1992| ) to indicate the presence of a circumstellar 
disk oriented perpendicular to the main outflow axis (Walther et 
al. 1993 , Torrelles et al. 1998), that presumably contains a massive 
protostar. 



1.2.3 Outline 

In this paper, we present the initial heterodyne data on Orion B 
taken for the GBS. We determine characteristic physical properties 
for two clouds, NGC 2024 and NGC 2071, comparing and con- 
trasting regions with varying physical characteristics. We present 
measurements of temperature and opacity across the regions, using 
these to provide temperature- and opacity-corrected masses and en- 
ergetics of the clouds and the high-velocity material. We use our 
13 CO data to investigate clumpy structure within the clouds, calcu- 
late the mass function of the condensations, and compare the values 
found with core mass functions in these regions, and to the initial 
mass function. Sec.plprovides an overview of the observations and 
data reduction. Sec.|3|describes the reduced data sets, highlighting 
regions of interest within the clouds with reference to previous ob- 
servations. Sec. [4] presents the physical characteristics of opacity, 
temperature, energetics and kinematics for the clouds, while Sec. [5] 
presents an analysis of the condensations as traced by 13 CO emis- 
sion. We summarize the results presented in the paper in Sec. [6] 



2 OBSERVATIONS 

2.1 Overview of HARP 

HARP is a recently-commissioned spectral-imaging receiver for 
the JCMT operating at submillimetre wavelengths ( [Buckle et al.| 
[2009l|Smith et al.|2003||2008] ). It works in conjunction with the new 
back-end correlator, ACSIS (Auto-Correlation Spectral Imaging 
System, [Buckle et al.|2009l|Dent et al.|2000| ) offering high-spectral 
and spatial resolution, the latter matched to that of SCUBA-2 ( |Hol-| 
|land et aL|2006| >. HARP maps spectral lines in the 325-375 GHz 
atmospheric window, the JCMT B-band frequencies, where tran- 
sitions of some of the most abundant molecules in the interstellar 
medium reside. The CO, 13 CO and C 18 / = 3 — > 2 lines in this 
window are used by the GBS to trace dense and/or warm gas around 
star-forming cores. The HARP imaging array consists of 16 SIS de- 
tectors arranged in a 4 x 4 grid, separated by 30 arc sec. The beam 
size is 14 arcsec at 345 GHz, which means this array arrangement 
under- samples the focal plane with respect to Nyquist and further 
data points must be taken in between the nominal on-sky positions 
to produce a fully-sampled map. HARP works in single sideband 
mode, with typical receiver noise temperatures under 150 K and 
main-beam efficiencies of T7 m b = 0.61 (Buckle et al. 2009). The K- 
mirror rotates the array with respect to the plane of the sky to max- 
imize the observing efficiency, for instance allowing large raster 
scan maps to be made with the array oriented at an angle to the scan 
direction. ACSIS provides either wide-band (up to 1.9 GHz wide) 
or high-resolution spectra (up to 31 kHz channels). Additionally 
the IF can be separated into two sub-bands to allow simultaneous 
imaging of two close-together transitions e.g. the 13 CO and C 18 
7 = 3^2 lines shown here. 

2.2 Description of Observations 

The observations presented here comprise 31.3 hrs of data taken on 
multiple nights from November 2007 to October 2008. The system 
temperatures, which are measured for each detector in every obser- 
vation, varied (across different weather bands) from 344 to 684 K 
for the 12 CO and 343 to 575 K for the 13 CO/C 18 data. Two fields 
in Orion B are analyzed: (i) NGC 2024, (10.8 x 22.5) arcmin 2 

centred on 05 h 41 m 39 s .8, -01° 54 2l" (J2000) at PA= -40 deg 
(east of north) (ii) NGC 2071, (13.5 x 21.6) arcmin 2 , centred on 

05 h 47 m 00 s .0, +00°19 / 53 // (J2000) angled at PA= -55 deg. The 
maps were taken in the raster position- switched observing mode, 
where spectra are taken "on-the-fly" with the telescope constantly 
scanning in a direction parallel to the sides of the map, taking spec- 
tra separated at 7.3 arcsec along this direction. The array is inclined 
at an angle of ~14 deg to the scan direction which produces rows 
separated by 7.3 arcsec perpendicular to the scan direction auto- 
matically ( Buckle et al. 2009). At the end of a scan row, the map is 
displaced half the array spacing perpendicular to the scan direction 
before a new row is started. This overlaps half of the new scan row's 
data points with the first to double the integration time and even 
out any noise variations due to missing receptors or intrinsic dif- 
ferences in receptor performance. The noise is further evened-out 
by observing a second map scanning in a perpendicular direction 
to the first to "basket-weave" the field. 

Separate absolute off positions were used for each field: 
(1) 05 h 39 m 00 s ,-01°00 / 00 // for NGC 2024 and (2) 05 h 43 m 44 s , 
+00°21 / 42 // .2 for NGC 2071. Each was verified to be absent of 
emission by examining short 60 s position- switched 'stare' ob- 
servations towards each position in 12 CO. A low-resolution 12 CO 
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Table 1. Current and required noise levels of GBS Orion B data 





current levels' 3 


required levels^ 


Source 


12 co 13 co c 18 co 


12 co 13 co c 18 co 




mean 1 a rms /K 


mean 1 a rms /K 


NGC 2024 


0.10 0.26 0.28 


0.3 0.25 0.3 


NGC 2071 


0.11 0.32 0.43 


0.3 0.25 0.3 



a 6 arcsec pixels, 16.6 arcsec beam, 1.0 ( 12 CO)/0.1( 13 CO/C 18 O) kms^ 1 
channels. 

b 1.5 arcsec pixels, 14 arcsec beam, 1.0 ( 12 CO)/0.1( 13 CO/C 18 O) kms^ 1 
channels. 



mode (1 GHz bandwidth with 977 MHz channel spacing) was used 
for this inspection as it provided lower noise and ensured the other 
weaker isotopologues would also be free of emission. 



2.2. 1 The GBS HARP/ACSIS configuration 

The 12 CO data were taken in an ACSIS dual sub-band mode: 

(1) with modest resolution to map the broad line flows, at a rest 
frequency of 345.796 GHz, with the 1 GHz bandwidth divided 
into 1024 channels, separated by 977 kHz (~0.85 km s _1 ) and 

(2) , simultaneously, with higher resolution to map the detail of 
outflows close to the ambient cloud velocity, again at a rest fre- 
quency of 345.796 GHz, with 250 MHz bandwidth divided into 
4096 channels, separated by 61 kHz (~0.05 km s _1 ). Similarly, 
the 13 CO/C 18 data were taken in a high-resolution dual sub- 
band mode, each sub-band has a rest frequency of 330.588 or 
329.331 GHz respectively with the 250 MHz band separated into 
4096 61 kHz channels (at a spacing of 0.055 km s- 1 ). 



2.3 Data Reduction 

The data were reduced using the Starlink project softwar^] with 
KAPPA (Currie et al. 2008 ) routines used to mask out poorly- 
performing detectors and calculate a self-flat (see Curtis, Richer & 
|Buckle|2009t for each observation. SMURF (Jenness et al.| 2008 ) 
routines were used to make the cube using a Gaussian gridding 
technique, with 6 arcsec pixels and a 9 arcsec beam, giving an 
effective angular resolution of 16.6 arcsec in the final datacubes. 
KAPPA routines were then used to remove a linear baseline and 
crop the edges of the cubes spatially and spectrally. These data are 
presented as a 'first look' at the data being gathered for the GBS 
towards Orion. Data in the isotopically-substituted species towards 
these regions are still being collected, and the final data sets will 
have a much increased signal-to-noise. Table [T] describes the mean 
noise achieved in the datacubes collected to date, and the noise re- 
quirements of the final GBS data products. Data are presented in 
units of antenna temperature (?2, |Kutner & Ulich |1981|), which can 
be converted to main beam temperature (T m b) using T m ^ = T£/ r/ m b. 



3 RESULTS 

3.1 HARP spectral imaging 

3.1.1 NGC 2024 

Fig. [2] displays the total integrated intensity images across the full 
width of the emission lines towards NGC 2024 in 12 CO 3^2 (- 
20.0-50.0 kms" 1 ), 13 CO 3^2 (4.0-16.0 kms" 1 ) and C 18 3^2 
(6.0-14.0 krns" 1 ). In this cloud, the integrated intensity image in 
12 CO is dominated by emission to the north of the map, particularly 
along a ridge extending south-west of the sources FIR 1-7 (marked 
in Fig. [2] with white crosses). A cavity is seen to the north-west of 
this position, with the rest of the 12 CO emission surrounding the 
cavity in a clumpy ring. The cavity marks a region of bright opti- 
cal emission. Fig. [3] shows the digitized sky survey (DSS) infrared 
imag^] overlaid with 12 CO contours. The dense ridge where we 
see the brightest molecular emission appears as a dark dust lane 
in the optical image. The bright optical filaments appear in regions 
which lack 12 CO emission. In particular, the large cavity contoured 
in 12 CO integrated intensity emission clearly follows the shape of 
the bright optical emission. The peak of the 13 CO emission is seen 
in the ridge containing the FIR sources. Further emission follows 
the 12 CO emission, delineating a cavity surrounded by a ring of 
clumpy, filamentary structure. The C 18 emission peaks on the po- 
sitions of the FIR sources in the ridge. Weaker, fragmented and 
clumpy emission surrounds the cavity seen in 12 CO and 13 CO. 

The 12 CO spectra show evidence of several velocity compo- 
nents along the line of sight. Fig. [4] shows the spectra from a sin- 
gle pixel at the peak of C 18 emission, from the three isotopo- 
logues. Also shown is a composite spectrum averaged over the ob- 
served region. The high-resolution 12 CO spectrum shows four dis- 
tinct peaks, at4.6kms _1 , 8.3 krns -1 , 11.2kms _1 and 13.0 kms -1 . 
The 13 CO spectrum has components at 4.2 krns -1 , 8.7 krns -1 
and 11.3 krns -1 , while the C 18 line has a main component at 
11.0 krns -1 , with a shoulder extension at 9.2 krns -1 . These line 
profiles suggest self absorption in the 12 CO line at 9.6 krns -1 and 
11.9 krns -1 . Previous observations have also found several ve- 
locity components towards NGC 2024 (e.g. Emprechtinger et al.| 
2009| . At 9 krns -1 , the dense ridge, thought to be in the fore- 
ground, is seen as a separate component in C 18 and 13 CO, while 
there is a dip in the 12 CO spectrum. This supports the suggestion 
that emission from the blue-side of the 12 CO line is being ab- 
sorbed by foreground materia l in the dense ridge (Emprechtinger 
et al.||2009||Graf et al.|l993| ). The main 13 CO and C ls O compo- 
nents peak near 11 krns -1 , the velocity of the extended molec- 
ular cloud in the background (e.g. |Barnes et al.||1989] >, and we 
also see a 12 CO component at this velocity. The dip in the 12 CO 
spectrum at 11.9 krns -1 , which is part of the main component of 
the 13 CO and C 18 emission, suggests that there may be a fur- 
ther foreground component at this velocity which is absorbing the 
12 CO emission. The self-absorption of the 12 CO emission undoubt- 
edly causes its different peak position to the isotopologues tracing 
higher column density material. The multiple line-of-sight compo- 
nents are not clearly seen in spectra averaged over the whole re- 
gion, and demonstrate the importance of high-spectral and spatial 
resolution observations in understanding the kinematics of star for- 
mation/molecular clouds, especially in clustered regions. 

Fig. [5] shows position- velocity (PV) diagrams of the three iso- 
topologues across a large section of NGC 2024. The integrated in- 
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from Gaia Skycat http://archive.eso.org/cms/tools-documentation/skycat 
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Figure 2. Integrated intensity images of 12 CO (top), 13 CO (middle) and 
C ls O (bottom) towards NGC 2024. Crosses mark the positions of sources 
FIR 1-7, and an isolated continuum source to the far south. Contours are 
12 CO: from 30 K kms" 1 with 20 K kms" 1 steps; 13 CO and C ls O: from 2 
K kms -1 , with contour levels doubling at every step. 
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Figure 3. DSS infrared image of NGC 2024 overlaid with 12 CO contours 
as in Fig. [2] 

tensity image in the left is marked with two lines showing the right 
ascension and declination cuts used. In the 12 CO map, we can see 
a velocity gradient running north-south, from ~13.0 kms -1 in the 
south to ~10.0 kms -1 in the north. High- velocity material is seen 
as extended filaments across much of the PV diagrams, with the 
brightest high- velocity material seen associated with the ridge in 
the declination cut. The cavity previously described is labelled in 
the right ascension cut, at two velocity intervals straddling the main 
bulk of the molecular cloud at 11.0 kms -1 . From 12.0-14.0 kms -1 
the cavity shows a distinct velocity gradient along its edge, and is 
narrower than the outline seen at velocities 8.0-10.0 kms -1 . The 
lower- velocity component of the cavity does not have the same ve- 
locity gradient along the edge. The very high- velocity outflow ma- 
terial is only clearly seen in the 12 CO PV diagrams, but emission 
from 13 CO shows extensions in velocity that suggest it is associ- 
ated with outflow material. Both the 12 CO and 13 CO show com- 
pact clumps of emission at blue-shifted velocities near 4 kms -1 , 
one of the components revealed in the spectra (Fig. [4]). The 12 CO 
PV diagrams show filaments extending between the clumps and the 
bulk of the emission, suggesting these clumps could be physically 
associated with the cloud, possibly through wind-driven activity. 
The emission at red-shifted velocities, from 12-14 kms -1 , shows 
a much smoother distribution than the blue-shifted emission. The 
C 18 emission is very compact in velocity, tracing the bulk of the 
dense material in the ridge, and no longer contains the component 
at 4 kms -1 . 

3.1.2 NGC 2071 

Fig. [6] shows total integrated intensity images across the full 
width of the emission lines for 12 CO 3^2 (-20.0-35.0 kms -1 ), 
13 CO 3^2 (-3.0-18.0 kms -1 ) and C 18 3^2 (5.0-12.0 kms -1 ) 
towards NGC 2071. The 12 CO emission is dominated by a curved 
filament of emission, surrounding a compact, bright, optical emis- 
sion region. Fig. [7] shows the DSS infrared image overlaid with 
12 CO contours towards NGC 2071, showing the 12 CO emission 
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Figure 5. Position velocity maps of NGC 2024 in 
image on the left. 



-CO, 



04:00 -2:00:00 56:00 52:00 -1:48:00 

Declination (ddd:mm:ss) 

13 CO and C ls O, through the positions shown in horizontal and vertical lines in the integrated intensity 



5:42:00.0 50.0 40.0 30.0 20.0 41:10.0 
Right ascension (hh:mm:ss.s) 



© 2009 RAS, MNRAS 000,[l]j20] 



GBS HARP Survey: Orion B 7 





Radio velocity (km/s) 

h oil! 

Figure 4. Spectra from a single pixel position (5 41 m 44 s .6 -1 55 33 , top) 
and averaged over the whole cloud (bottom) of 12 CO (solid line) 13 CO 
(dashed line) and C ls O (dotted line) towards the C ls O peak in NGC 2024. 
All of the spectra have been observed at high-resolution, and binned to 
0.1 kms" 1 . 



nearly surrounding the bright nebula. The 13 CO and C 18 emission 
peaks along the compact ridge, but both show a different spatial 
structure to the 12 CO emission at the southern end of the ridge, at 
the position of the optical nebula. In 13 CO and C 18 0, the emission 
to the south, where the 12 CO integrated intensity falls off, displays 
separate arcs. 

The spectra towards individual positions in NGC 2071 are 
shown in Fig. M At the systemic velocity of the cloud, around 
10 kms -1 , the CO is heavily self-absorbed, even in the low- 
density regions near the optical nebula (Fig. [8] top). Towards the 
peak of emission from C 18 (Fig.[i] middle), emission from 13 CO 
is also self-absorbed. The absorption dip in 12 CO emission is also 
clear in the spectra averaged across the whole cloud (Fig. [8] bot- 
tom). 

Fig.|9]shows the PV diagrams for NGC 2071, as for Fig.|5] We 
see high- velocity material from 12 CO and 13 CO, and evidence for 
the C 18 emission tracing lower- velocity outflow material. There 
is no clear velocity gradient through the cloud. The 12 CO and 13 CO 
PV diagrams along a cut in right ascension suggest that there may 
be a second component along the line of sight to the west, or that 
the cloud is more extended in velocity in this region. In contrast to 
NGC 2024, towards NGC 2071, we do not see any clear features 
m the 12 CO PV dia gram that may be associated with cavities or 
bubbles carving out to the cloud edges through activity within the 
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Figure 6. Integrated intensity images of 12 CO (top), 13 CO (middle) and 
C ls O (bottom) towards NGC 2071. Contours are 12 CO: 30,35 (white) and 
45,65,105,185,345 (black) K kms" 1 ; 13 CO 7,9 (white) and 11,15,21,29,39 
(black) K kms" 1 ; C ls O: 2 (white) and 4,7 (black) K kms" 1 . The first C ls O 
contour is at 5a. 
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Figure 9. Position velocity maps of NGC 2071 in 12 CO, 13 CO and C 18 Q, as for Fig.[5] 





cloud. The C 18 emission shows a compact, clumpy structure in 
the PV diagrams. 



3.2 High- velocity material 

Figs.[To|and[TT|show false colour red, green and blue images of the 
average intensity across the 12 CO line profile towards NGC 2024 
and NGC 2071, where the emission near the systemic velocity of 



the clouds is shown in green, and the red- and blue-shifted gas is 
shown in the corresponding colour. 

Towards NGC 2024, a single region dominates the high- 
velocity red-shifted material, associated with an extended outflow 
from FIR5 and a compact outflow from FIR6 (Richer et al. 1989 ). 
Additionally, more diffuse components are seen throughout the 
map, which probably indicate different components along the line- 
of- sight, rather than protostellar outflows. However, there are two 
regions to the north, and several regions in the centre of the map 
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Figure 7. DSS infrared image of NGC 2071 overlaid with 12 CO contours 
as in Fig. [6] 

where the relative position of red and blue emission, and the shape 
of the emission, suggest there are embedded outflow sources. The 
PV diagrams (Fig. [5} show several filamentary structures in veloc- 
ity that may be associated with outflows that are less extensive in 
velocity than the obvious extremely high- velocity flows. 

Fig.[T2]shows contours of blue- and red-shifted emission over- 
laid on the integrated C 18 intensity map for the main NGC 2024 
outflow source. In the 12 CO J = 3 emission, we see the high- 
velocity blue- shifted lobe is spatially compact, although extended 
in velocity. This is the only region in NGC 2024 where such high- 
velocity emission is detected. 

NGC 2071 is dominated by a high- velocity outflow, previ- 
ously identified as one of the most energetic and collimated out- 
flows ever discovered (e.g. Stojimirovic et al. 2008). The false 
colour image shows there are at least two more outflows, as sug- 
gested by Stojimirovic et al. ( 2008). The PV diagrams (Fig. [9} sug- 
gest how powerful this outflow is, with small extensions in velocity 
showing in the C 18 emission. Again, the velocity-extended fila- 
mentary structures in this large-scale PV diagram suggest that there 
may be several protostellar outflows embedded within the cloud. In 
Sec. |4.5[ we investigate the energy injection from outflows in these 
regions. Fig[T3|shows the blue- and red- shifted contours overlaid on 
theC 18 emission for the main NGC 2071 outflow region, overlaid 
with SCUBA sources (J06). Red- and blue-shifted outflow lobes 
can be identified, for at least three, and possibly five sources. 

3.3 Low density material 

Towards both regions, a cavity is seen in emission from 12 CO, asso- 
ciated with an optically-bright Hll region. However the behaviour 
of the less abundant isotopologues, and the velocity structure of 
these regions is different towards the two clouds. 

Towards NGC 2024, emission from all of the CO isotopo- 
logues falls off in this region, as can be seen in the integrated inten- 
sity maps (Fig. [2]). The cavity is coincident with the bright optical 
emission, and it may be that we are seeing the CO disassociated by 
an Hll region. F ig.[T4lshows the SCUBA 850 \xm continuum image 
( |Di Francesco et alj 2008 ) overlaid with C 18 integrated intensity 
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Figure 8. Spectra from single pixel positions; (top) 5M7 m 04 s .4 +0 14 53 

, o I II 

towards the optical nebula; (middle) 5M7 m 04 s .O +0 21 53 , towards the 
C ls O peak; and averaged over the whole cloud (bottom) of 12 CO (solid 
line) 13 CO (dashed line) and C ls O (dotted line) in NGC 2071. 



contours. The C 8 emission closely follows the structure seen in 
the dust, and the dumpiness of both the dust and the C 18 emission 
suggest that cores are present in a region surrounding the cavity. We 
present an analysis of all the 13 CO condensations in Sec. [5] 

In the NGC 2024 12 CO PV diagram (Fig. [5), the cavity can 
clearly be seen, and shows a different spatial structure at high and 
low velocities. There is an abrupt cut-off in emission across the 
lower velocities, from 8.0 to 10.0 kins" 1 . At the higher velocities, 
from 12.0 to 14.0 kms -1 , there is a velocity gradient, and the east- 
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Figure 10. False colour image of the average 12 CO emission towards NGC 
2024 in red (13.2-19.2 kms -1 ), green (7.2-13.2 kms -1 ) and blue (1.2- 
7.2 kms" 1 ) velocity channels. 
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Figure 11. False colour image of the average CO emission towards NGC 
2071 in red (13.1-19.1 kms -1 ), green (7.1-13.1 kms" 1 ) and blue (1.1- 
7.1 kms -1 ) velocity channels. 
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Figure 12. NGC 2024 high-velocity red- and blue-shifted 12 CO emission. 
Red contours show emission from 14.5 to 40.0 kms -1 , and blue contours 
show emission from -20.0 to 4.0 kms -1 , both starting at 5 Kkms -1 , in 
steps of 20 Kkms -1 . The image shows the integrated C ls O intensity. The 
suspected driving sources FIR5 and FIR6 have been labelled. 
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Figure 13. NGC 2071 high-velocity red- and blue-shifted 12 CO emis- 



-l 



blue 



sion. Red contours, show emission from 14.5 to 35.0 kms 
contours show emission from -25.0 to 0.5 kms -1 , contour levels at 
5,10,15,25,45,65,85,105,125,145 Kkms -1 . The image shows the inte- 
grated C ls O intensity. Black circles denote positions of SCUBA cores (J06). 



ern edge has a higher gradient in velocity than the western edge. 
If we are seeing an Hll region, then the different morphologies 
between the gas at low velocities and at high velocities may be 
giving information on the three-dimensional spatial structure. The 
low velocities, tracing material at the front of the cloud, have an 
abrupt cut-off to the emission because at this point, the material 
has broken out of the molecular cloud, and there is no more mate- 
rial entrained in the boundary layers. At higher velocities, if tracing 
material deeper into the cloud, the eastern edge is pushing into the 
dense molecular ridge seen in the C 18 emission, while the west- 



ern edge is pushing into more clumpy and fragmented material. 
This would account for the different velocity gradients seen in the 
PV diagrams, and supports the three-dimensional models suggested 
for this region ( Emprechtinger et al. 2009, Barn es et al.|1989| >. 

Fig. [15] shows the NGC 2071 SCUBA 850 |xm map <[Di] 



2008), overlaid with C 18 integrated intensity con- 



Francesco et al. \ 

tours. The C lb O emission does not follow the structure seen in the 
dust as closely as in NGC 2024, and it is in the region co-incident 
with the optical nebula that we see the most discrepancy. In the 
PV diagrams (Fig. [9]), we see that the 13 CO and C ls O emission is 
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Figure 14. NGC 2024 SCUBA 850 jam maps (Di Francesco et al. 2008 1, 
overlaid with C ls O integrated intensity contours. Contour levels are as 
Fig.0 
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Figure 15. NGC 2071 SCUBA 850 f^m map \Di Francesco et al.| 2008 1, 
overlaid with C ls O integrated intensity contours. Contour levels are as in 
Fig.© 



clumpy and fragmented, both spatially, and in velocity, while the 
12 CO emission has a lower intensity region near the systemic cloud 
velocity, embedded within material at much higher velocities. 



4 ANALYSIS OF CLOUD PARAMETERS 
4.1 Opacity 

We calculate the gas opacity based on the line peak ratios of 
12 CO/ 13 CO and 13 CO/C 18 0. The intensity ratio is related to the 
opacity through: 



( 1318 



r;( 13 co) _ i-ex P (-T 13 ) 



r;(ci8 ) i-ex P (-T 18 ) 



(i) 



= r;( 12 co) i-ex P (-T 12 ) 

1213 r;(C 13 ) l-exp(-T 13 ) 



(2) 



This method makes several assumptions about emission from 
the isotopologues (see, for example Myers, Linke & Benson 1983 ; 
|Ladd, Fuller & D eane 1998 for details). In particular that the beam 
efficiency, filling factor, and excitation temperature are the same for 
all the isotopologues, and that the ratio of line widths to the line of 
sight extent of the emitting gas is the same for all isotopologues. 

Where the lines are optically thin, the ratio R should tend to 
the abundance ratios, X[ 13 CO/C 18 0]=8.4 an d X[ 12 CO/ 13 CO1=70 
for 13 CO/C 18 Q and 12 CO/ B CO respectively {Frerking, Langer & 
|Wilson| 1982] |Wilson|l999l >. The 12 CO/ 13 CO intensity ratio maps 
for each source are shown in Fig. [T6] 

From the ratio maps, we can estimate the abundance ra- 
tios from the values the ratios tend to towards the edge of de- 
tected emission regions, where we expect the lines to be optically 
thin. At the edges of the C 18 emission regions, Rms tends to 
8-12 in NGC 2024, and we therefore adopt X[ 13 CO/C 18 O]=10 
for NGC 2024. For NGC 2071, R ms tends to 6-9, so we 
similarly adopt X[ 13 CO/C 18 0]=7.5. We detect 12 CO and 13 CO 
across the full extent of the mapped region, and the integrated 
intensity ratio does not approach the expected abundance ratio 
of —70. The highest ratios that we detect are —20, suggest- 
ing that 12 CO is optically thick in the line core everywhere in 
both clouds, and so we use X[ 12 CO/ 13 CO]=70. Using the above 
equations and abundance ratios, mean opacities in the clouds are 
f( 12 CO) -119, f( 13 CO) -1.7, and f(C 18 0) -0.17 in NGC 2024, 



while f( iZ CO) -157, f( ij CO) -2.24, and f(C i8 0) -0.3 in NGC 
2071. These values suggest that C 18 is generally optically thin, 
while 13 CO is marginally optically thick, and 12 CO is optically 
thick throughout both clouds. Note that self-absorption in the line 
profiles, which we see towards the densest regions, means that the 
opacity will be over-estimated for 12 CO, and to a lesser extent for 
13 CO. 



4.2 Temperature 

Assuming local thermodynamic equilibrium (LTE), the excitation 
temperature (T ex ) of the gas can be calculated from line peak tem- 
perature (r max ) using the standard radiative transfer relation in an 
isothermal slab. Also assuming that the 12 CO emission is optically 
thick (t — ► 00) and fills the beam in the line core across the cloud, 
then (e.g. |Pineda, Caselli & Goodman|2008] >: 



r eX (3^2): 



16.59K 



In [1 + 16.59K/ (r max ( 12 CO) + 0.036K)] ' 



(3) 



for regions where the 12 CO line is not self-absorbed. We can iden- 
tify the regions where there is self- absorption by comparison of 
the 12 CO line profile with that of the 13 CO line profile. Where the 
12 CO is optically thick and self-absorbed, the 13 CO line peak tem- 
perature is higher than the 12 CO line peak temperature. Towards 
these regions, 13 CO is optically thick, and we use the equivalent 
form of Eq.[3]to calculate T QX from 13 CO. Fig. 17 shows Tq X across 



NGC 2024 (top) and NGC 2071 (bottom), calculated using both 
12 CO and 13 CO. 

NGC 2024 generally reaches higher excitation temperatures 
than NGC 2071. Towards NGC 2024, the mean excitation temper- 
ature is 31.8 K, which increases to >50 K in the dense regions, and 
—70 K in regions which are known to contain energetic outflows. 
Towards NGC 2071, the mean excitation temperature is 19.6 K, 
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Figure 16. The 12 CO/ 13 CO (left) and 13 CO/C 18 (right) peak intensity ratios for regions with detections <3a towards NGC 2024 (top) and NGC 2071 
(bottom), overlaid with C ls O integrated intensity contours. The 13 CO/C 18 ratio data has been smoothed by a Gaussian with 12 arcsec FWHM to make the 
features more visible in these images. 



which increases to >30 K in the dense regions, and ~50 K in the 
regions associated with energetic outflows. 

Towards the dense cores of Orion B, kinetic temperatures have 
been published using line ratios of H2CO, a molecule known to 
trace gas kinetic temperatures ext remely well {[Mangum & Woot-| 
|ten|1993||Tothill & Mitchell|2001||Mangum et al.|1999| ). Although 
only a few of these sources overlap with the regions we have ob- 
served, the excitation temperatures calculated from CO agree well 
with those calculated from H 2 CO, except towards NGC 2024 FIR5, 
the probable driving source of the energetic outflow seen as a large 
red-shifted outflow lobe (Fig. [12) . Average values for the small 
number of H2CO dense cores, with the exception of NGC 2024 
FIR5, are 44 K from H 2 CO, and 39 K from CO. Given that the 
H2CO transitions have a critical density of ^3.9 x 10 6 cm -3 , while 
the critical density of the 12 CO transition is ~3.5x 10 4 cm -3 , and 
that the velocity structure of these regions is relatively complex, 
this is surprising, and suggests the combination of 12 CO and 13 CO 
is a relatively good tracer of excitation conditions in dense re- 
gions. We have also compared our CO excitation temperatures with 
those derived from dust temperatures towards SCUBA cores in both 
clouds by J01, J06. Fig. [T8] shows the integrated intensity maps of 



CO, overlaid with contours of C O, and marked with the posi- 
tions of the SCUBA dust cores from J01, J06. The cores are marked 
in green if T QX ( l2 CO)/Td ust <1.5, and in red if the ratio is >1.5. For 
the small number of cores that also have H2CO measurements, the 
average dust temperature is 41 K, similar to the H2CO and CO val- 
ues. Towards NGC 2071, all of the SCUBA cores have CO and 
dust temperatures ratios <1.5, while towards NGC 2024, the ra- 
tio is <1.5 only towards the high column density regions of NGC 
2024, along the molecular ridge. Towards NGC 2024, the cores 
where dust temperatures are similar to temperatures derived from 
CO are spatially grouped towards the highest column density re- 
gions. In NGC 2024, the differences in r ex ( 12 CO) and T^ ust towards 
the lower column density regions may be evidence of sub-thermal 
excitation of CO. Given the complex velocity structure of these re- 
gions, and the complications that opacity adds to the line profile, 
observations utilizing all CO isotopologues so that opacity correc- 
tions can be made appear to be relatively good tracers of the gas 
and dust temperatures in energetic, densely clustered regions. 
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Figure 17. Excitation temperature (in K) for NGC 2024 (top) and NGC 
2071 (bottom). To indicate regions where the 12 CO may be self-absorbed, 
and so underestimating the excitation temperature, the contours show C ls O 
integrated intensity. 



4.3 Mass and energetics 

Since we have isotopic data, we can calculate the mass in the large- 
scale cloud, using the 12 CO and 13 CO data, which we detect ev- 
erywhere. In the dense molecular ridge, where 12 CO and 13 CO line 
profiles are affected by self-absorption, we can use the C 18 data, 
assuming LTE conditions. We compute the masses assuming T ex is 
the same for all three isotopologues, and equal to the mean 12 CO 
excitation temperature (Sec. |4.2| ). We correct for optical depth ef- 
fects in the 12 CO a nd 13 CO data usi ng the mean opacitie s f ( 13 CO), 
and f( 12 CO) (Sec.|4~I}. Following |Garden et al.|fl99l} : 



N(C ls O 



N( 13 CO) 



8.26 x 10 13 exp 
r ex + 0.88 



15.81 



1 — exp( : 



88 r 



8.28 x 10 13 exp 



rdv cm 
15.87 N 



7k + 0.88 r 
-exDf^W 



1 — exp(— £ 



rdv cm 



(4) 
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Figure 18. 12 CO integrated intensity maps for NGC 2024 (top) and NGC 
2071 (bottom), overlaid with C ls O contours. SCUBA cores from J01, J06 
are marked in green where the ratio r e x( 12 CO)/7d us t ^1.5 and in red other- 
wise. 



N( 12 CO) 



7.96 x 10 13 exp( — 
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Tdv cm 
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with Tqq = 2.7K. The mass is given by: 



Mi 



gas • 



1.13 x 10 - V// 2 m#d 2 AaA/3XAf(CO) 



(11) 



where AaAfi is the pixel area in arcsec 2 , d is the distance in 
pc, v is in km s -1 and X is the isotopic abundance ratio of CO 
and its isotopologues relative to H2. We adopt a mean molecu- 
lar weight per H2 molecule of \ih 2 = 2.72 to include helium, and 
X[ 12 CO]=1.0x 1(T 4 , X[ 13 CO]=1 .4x 1(T 6 , X[C 18 0]=1.7x 1(T 7 
( |Frerking et al.|1982[|Wilson|19991 ). The results are listed in Tab.[2] 
The masses that we find for optically thick 12 CO are smaller than 
previous mapping observations of a 19 deg 2 region in Orion B, us- 
ing a lower excitation transition, where a total mass of 8.3 x 10 4 M 
was measured ( Maddalena et al. 1986 ). The total mass we derive for 
NGC 2024 and NGC 2071, covering a region -100 times smaller 
in area and using optically thick 12 CO, is 8.5 xlO 3 M . In 13 CO, 
we find a total mass of 2.5 xlO 3 M . The total mass of gas that 
we detect in C 18 towards both clouds is 1000 M . Emission from 
the optically thin isotopologue C 18 traces a smaller region within 
the cloud, and does not suffer from the same uncertainties due to 
optical depth corrections. If the opacity corrections using mean val- 
ues were accurate, we would perhaps expect the masses calculated 
from 12 CO and 13 CO emission to be similar to that calculated from 
C 18 emission. However, the differences may be due to physical 
conditions within the clouds affecting the transitions differently. 
Towards NGC 2024 in particular, the many velocity components 
seen in the spectra indicate that a single excitation temperature and 
density along the line of sight may not be a good approximation of 
the conditions. The clouds may not be in LTE, and emission from 
the isotopologues subthermally excited. As suggested by Pineda 
|et al.| ( |2008| ), the shielding of material within the cloud leads to 
physical and chemical properties which can affect 12 CO and 13 CO 
differently to the rarer isotopologue C 18 0. We find a total mass 
of 1111 M using 13 CO data, if we assume that 13 CO emission 
is also optically thin. The total mass calculated from 12 CO under 
the assumption of optically thin emission is 65 M , indicating that 
emission from this isotopologue is, as calculated above, very opti- 
cally thick, and tracing only the surface regions of the clouds. NGC 
2024 contains more mass than NGC 2071, measured in all of the 
isotopologues. 

The gravitational potential energy, W, of the cloud can be 
calculated assuming a density distribution p r a (M acLaren, 
Richardson & Wolfendalej 19881 |Williams, de Geus & Blitz| 19941 



Bertoldi & M cKee 1992). The kinetic energy, E\^ n , can be cal- 
culated by estimating the three-dimensional velocity dispersion, 
°v 3D' f rom me one-dimensional velocity dispersion of the average 
spectrum, Cu co , 0"i 3co or Oc 18 o : 



W 
1 



3 GMi 



5' R 

l+a/3 
l + 2a/5 
1 ? 



J v3D 



J CO~ 



kT 
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m C o 



(12) 
(13) 
(14) 
(15) 



where G is the gravitational constant, k is Boltzmann's constant, 
T is the kinetic temperature of the gas, which we have taken to 
be equal to the mean 12 CO excitation temperature and mco is the 
atomic mass of 12 CO, 13 CO, or C 18 0. We assume the clouds follow 
p oc r -2 density distributions, yielding 7=5/3. The LTE masses, 



kinetic and gravitational potential energies of both clouds are listed 
inTab.[2] 

NGC 2024, as well as containing more mass than NGC 2071, 
has greater kinetic and gravitational potential energies. The gravi- 
tational potential energy calculated from 12 CO and 13 CO is larger 
than the kinetic energy in both regions, implying that these clouds 
are both gravitationally bound. The emission from C ls O, tracing 
much smaller regions containing all of the currently active out- 
flows, has more kinetic energy. The line profiles and PV diagrams 
of 12 CO and 13 CO indicate that both of these isotopologues are 
tracing some of the high-velocity, presumably outflowing, mate- 
rial. In Sec. |4.5[ we compare these values to the energetics in the 
high- velocity material. 

4.4 Cloud kinematics 

Fig.[T9]shows the velocity at the line peak of the optically-thin gas. 
This has been created from the 13 CO and C 18 data, using 13 CO 
in regions where the C 18 line peak falls below the 3(7 level. To- 
wards NGC 2024, a velocity gradient runs across the image, from 
8 kms -1 in the north, to 10 kms -1 in the south, with the molecular 
ridge in the south seen in emission at 11.5 kms -1 . The dust cavity 
is seen in material at > 1 1 .5 kms -1 , surrounded by material at 8.0- 
10.0 kms -1 , with a further higher velocity band to the northwest. 
In the southernmost part of the map, a band of high- velocity mate- 
rial appears, which is the tip of the NGC 2023 region. The velocity 
gradient can also be clearly seen in the 12 CO PV diagram (Fig.|5J). 
Towards NGC 2071, there is no obvious underlying velocity gradi- 
ent. All of the velocity variations seen in the NGC 2071 map appear 
to be associated with current star forming activity in the cloud. In 
the south, there is a large region that is clearly red- shifted with re- 
spect to the rest of the cloud, which is in the same direction as the 
very energetic red-shifted outflow lobes seen towards NGC 2071 
(Fig.[l3j 



Fig. 



20 



) shows the equivalent width (/ Tdv/Tp^) for the opti- 
cally thin emission towards the two regions. The exterior edges of 
the molecular ridge towards NGC 2024, traced by the outer con- 
tours of C 18 0, have lower velocity dispersions than the material 
within which it is embedded, and the energetic region near the FIR 
sources. Towards NGC 2071, the material with a lower velocity 
dispersion is perpendicular to the outflow direction. In the regions 
extending along the outflow, a cone- shape of higher velocity dis- 
persion material can be seen. 



4.5 Mass and energetics of the high- velocity material 

Using the line centre of the average C 18 spectrum to identify 
the velocity of the cloud, we calculate the total mass in the high- 
velocity material generated from the outflows using 12 CO emis- 
sion, with red- and blue-velocity extents across the line profile. 
We assume that the high-velocity emission is optically thin, and 
at a similar excitation temperature to the bulk of the gas, given 
by the mean 12 CO excitation temperature, following the method 
in Sec. |4.3| We obtain masses of a few x 1 M . This is much 
lower than found using lower excitation transitions of 12 CO and 

although 



13 CO (e.g. |Snell et al.||l984[ |Stojimirovic et al. 



2008 



these authors used a combination of isotopologues across the ve- 
locity range for their calculation). For 13 CO, where the assumption 
that the emission is optically thin in the line wings is more likely 
to be correct, we obtain masses ~ a few x 10 M , in better agree- 
ment with those previously found, which suggests the optically thin 
assumption for 12 CO line wings produced in outflows is not valid. 
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Table 2. Table of masses and energetics 



Cloud 


LTE mass(xlO 3 /M ) 


£ kin (xl0 40 /J) 






W(xl0 40 /J) 






12 co 13 co c 18 o 


12 co 13 co c 


18q 


12 co 


13 co c 


18q 


NGC 2024 


4.8 1.6 0.6 


37.5 5.9 


1.6 


9.5 


1.1 


0.2 


NGC 2071 


3.7 0.9 0.4 


24.5 1.6 


0.4 


5.7 


0.3 


0.1 




50 40 30 20 10 5:47:00 50 40 30 20 46: 1( 
Right ascension 

Figure 19. NGC 2024 map of the velocity at the line peak (in kms -1 ) for 
C ls O where emission is above the 3a level, and 13 CO at all other regions 
(top). 13 CO is detected above the 3a level at all points in the map. The same 
Fig. for NGC 2071 (bottom). 



From the 3 CO emission, we calculate the energy and momen- 
tum of this high velocity material using : 

P = MV char (16) 
Ek = \mvI w (17) 

where V c h a r is a characteristic velocity estimated as the difference 
between the maximum velocity with detectable emission in 12 CO 
emission and the cloud velocity identified from the line centre ve- 
locity of a Gaussian fit to the average C 18 line profile. For this 
reason, and since we do not account for the inclination to the line of 
sight for any outflows, the calculated values are upper limits. The 




20 10 5:42:00 50 40 30 20 10 41:00 



Right ascension 
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Right ascension 

Figure 20. NGC 2024 map of equivalent width (in kms -1 ) for C ls O where 
emission is above the 3a level, and 13 CO at all other regions (top). 13 CO 
is detected above the 3a level at all points in the map. The same Fig. for 
NGC 2071 (bottom). Data has been smoothed by a Gaussian with 12 arcsec 
FWHM to make the features more visible in these images. 



mass, momentum and energy values for the red- and blue-shifted 
emission for both clouds are listed in Tab. [3] 

The total energy of the high- velocity material within the cloud 
is 100 times more than that calculated by investigating individual 
flows. Towards NGC 2071, |Stojimirovic et al.| ( [2008] ) found a total 
energy in the red and blue outflow lobes of 9.67 x 10 45 ergs, com- 
pared to our value of 1 1.15 x 10 47 ergs. The values we calculate are 
likely to be over-estimated, since we have used a maximum esti- 
mate for the characteristic velocity. Comparing these values to the 
kinetic energy of the cloud as measured from 13 CO emission, the 
high-velocity material accounts for very little of the total kinetic 
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Table 3. Mass and energetics of high-velocity material 



Cloud Velocity 


Vchar 


CO Mass 


J CU Mass 


bnergy 


Momentum 




kms" 1 


M 


M 


xl0 40 J 


M kms -1 


NGC 2024 red 


32.3 


4.50 


19.64 


2.03 


634 


NGC2024 blue 


25.2 


3.38 


14.59 


0.92 


367 


NGC 2071 red 


43.0 


6.75 


37.34 


6.86 


1604 


NGC 2071 blue 


39.2 


5.66 


28.01 


4.29 


1099 



energy in NGC 2024, while the high-velocity material dominates 
the kinetic energy in NGC 207 1 . 



5 ANALYSIS OF THE GASEOUS CONDENSATIONS 

The structure of molecular clouds has been characterized in numer- 
ous ways in the literature, including identifying discrete objects in 
the emission (often termed clumps). By studying the properties of 
such clumps we hope to explore the properties of the cores, har- 
boured by some of the clumps, which will go on to form stars. Fre- 
quently, automated clump-finding routines have been used, to look 
for clumps in an unbiased manner, in both spectral-line and dust 
continuum datasets. In this section, we undertake a statistical de- 
composition of the structure of the 13 CO 7 = 3^2 emission using 
the CLUMPFIND algorithm ( Willia ms" et al.|1994), implemented in 
Starlink's new CUPID package (Be rry et al.|2 007 ). The optically- 
thin C 18 0/ = 3^2 data will be a better tracer of the mass in 
the star-forming cores, where the 13 CO line is saturated. Here we 
explore the 13 CO 7 = 3^2 bulk cloud emission. Many different 
terminologies have been used to name clumps in molecular clouds, 
here we call any individual object identified using CLUMPFIND a 
'condensation', to emphasize that it may not collapse to form a 
star. 

The data were smoothed to 0.25 kms -1 spectrally to improve 
the signal to noise ratio but still maintain sufficient resolution to 
identify condensations with the narrowest anticipated line widths. 
At our assumed cloud temperatures (50 and 30 K for NGC 2024 
and NGC 2071 respectively) the 13 CO line widths expected from 
thermal broadening alone (y/kT jmnmco) are 0.4 kms -1 for NGC 
2024 and 0.3 kms -1 for NGC 2071. At this resolution, the mean 
RMS noise across the map is reduced to 0.14 and 0.20 K for NGC 
2024 and NGC 2071 respectively. We ran CLUMPFIND using a 
level spacing of 2orms (the value recommended by Willi ams~et| 
|aT]([r994 ) to reduce the contamination by noise) and a lowest level 
of 2orms> which ensures the minimum peak intensity of a con- 
densation is 4orms- Additionally, condensations were rejected if 
they touched any edge of the data array, contained fewer than 16 
(3D) pixels or were smaller than the beam size. These parameters 
resulted in a catalogue of 1561 and 1399 condensations in NGC 
2024 and NGC 2071 respectively. 

Previous studies of the clump distribution in the dust contin- 
uum emission of Orion B have identified a factor of ~ 10 fewer 
condensations than we do from the 13 CO emission (J01, J06). This 
is probably because the SCUBA observations are not sensitive to 
flux on large-scales whereas the 13 CO gas should trace this more 
ambient emission, which can account for the 13 CO condensations 
identified which are not spatially associated with a dust clump. Ad- 
ditionally, the spectral information in the 13 CO data may break up 
multiple condensations at distinct velocities along the line of sight, 
which may be superposed in continuum data. For each dust clump 
identified by J01 and J06 in our observing fields, we found the mean 



number of 13 CO condensations whose peak lies within the radii of 
the dust clumps to be - 7.4 for NGC 2024 and - 3.7 for NGC 
2071. 



5.1 Condensation properties 

In the catalogue generated by CUPID, the 'size' of each condensa- 
tion along each axis (right ascension, declination, velocity) is given 
by the standard deviation of the pixel co-ordinate values about the 
centroid position, weighted by the pixel values, then corrected to 
remove the effect of instrumental smoothing. We define the radius 
of a condensation as the geometric mean of the size along axes 1 
and 2(fli andtf 2 ): 

R=^R 2 , (18) 

where R\ and R 2 have been already been deconvolved with the 
beam size. 

For NGC 2024 the condensation radii range from 1.3 x 
10 -3 pc to 0.053 pc with a mean value of 0.018 ±0.009 pc. For 
NGC 2071 the values are similar; the radii range from 1 .3 x 10 -3 pc 
to 0.075 pc with a mean value of 0.017 ± 0.008 pc. 

We estimate the three-dimensional velocity dispersion of each 
condensation from the one-dimensional 13 CO velocity dispersion, 
<7i3 CO (the size along axis 3 calculated by CUPID), as described in 
Sec. |4.3| For NGC 2024, the velocity dispersions range from 0.66 
to 2.0 kms -1 , with a mean value of 0.99±0.23 kms -1 . Again the 
values are similar for NGC 2071, ranging from 0.50 to 2.5 kms -1 
with a mean of 0.88±0.19kms -1 . 

Assuming the clouds are in LTE and the 13 CO 7 = 3^2 tran- 
sition is optically thin, the masses of the 13 CO clumps can be de- 
rived following Sec. |4.3| At moderately high densities (which we 
expect the 13 CO / = 3 — > 2 transition to probe, given its critical 
density) the cloud gas and dust should be coupled and at the same 
temperature (Burke & Hollenbach 1983 ). The mean temperatures 
of the dust clumps associated with our 13 CO condensations, com- 
puted by J06 and J01, are 45 ± 1 1 K and 26 ± 8 K for NGC 2024 
and NGC 2071 respectively. We therefore adopt excitations tem- 
peratures of T ex = 50 and 30 K for NGC 2024 and NGC 2071 re- 
spectively in our mass calculations, which are also consistent with 
the temperatures we derived in Sec. |4.2| (see Fig.fTT]). 

In NGC 2024, we find the masses of the condensations 
range from 2.6 x 10 -3 M to 13 M©, with a mean value of 
0.52 ± 1.2 Mq. For NGC 2071 the condensations are less mas- 
sive, ranging from 3.4 x 10 -3 M to 5.3 Mq, with a mean value 
of 0.22 ± 0.34 M . These masses are much less that those derived 
from the dust analysis of J06 and J01, who found clumps as massive 
as 90 M in NGC 2024 and 30 M in NGC 2071. As mentioned 
above, this could be because there are several 13 CO condensations 
associated with each dust clump. Furthermore, opacity effects and 
freeze-out of CO in dense regions could reduce the observed inten- 
sity and hence the mass estimate. 
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The virial mass of a spherical condensation with a density pro- 
file of p oc r a is calculated using ( [MacLaren et al.|1988[|Wilriams| 
|etal.|1994) : 



from the regressio n coefficient in log-log space. Ikeda, Kita muraj 



37G 



(19) 



where R is the radius of the condensation, o v ^d is the three di- 
mensional velocity dispersion (Eq[T5]> and G is the gravitational 
constant. Assuming the condensations have density distributions of 
p oc r -2 , 7=5/3. The virial masses for NGC 2024 range from 0.21 
to 30 M with a mean of 5. 1 ± 4.6 M©, while for NGC 2071 M vir 
ranges from 0.13 to 32 M with a mean, 3.4 ±2.7 M . The virial 
masses are much higher than the LTE masses for every condensa- 
tion, but this is partly because the LTE masses are underestimates. 



5.2 Correlation between the condensation properties 

In Fig. [21] we plot the Mlte — R and ov — R relations for both 
clouds. In the plots of Mlte versus R, there appears to be a sep- 
arate population of low-mass condensations. These low-mass con- 
densations are roughly separated by the line Mlte/M© = R/pc and 
are plotted in red. The remaining condensations (plotted in black) 
appear to have a strong correlation between mass and radius, with 
correlation coefficients of 0.8 for both clouds. The linear regression 
coefficient (s xy /s xx ) of log (Mlte/M©) versus \og(R/pc) is 2.6 for 
NGC 2024 and 1.7 for NGC 2071, implying that the mass-radius re- 
lationship is of the form Mlte/M© oc (fl/ pc )~ 2 - 6 and oc (R/pc)~ lJ 
for NGC 2024 and NGC 2071 respectively. Given the scatter in the 
plots and the uncertainties in the values of Mlte an d R, the rela- 
tion for both clouds is consist ent with the one found by |Kramer,| 
Stutzki & WinnewisserlJ 19961 for their 13 CO condensations in the 



southern part of Orion B, which had a power law index of 2.2. The 
relations are also consistent with Larson 's Law relatin g mass and 



21 



radius, which is of the form Mlte 00 

R 2 ( |Larson|l98l] ). In Fig. 
the linear regression fits are shown with the solid black line, and 
the Larson relation is shown with the broken line. 

The distinct population of low-mass condensations, on the 
other hand, only have a weak correlation between mass and radius 
(with correlation coefficients of 0.5 in both clouds). Since these 
populations have such low masses and behave differently to the 
rest of the condensations it is possible they have some contamina- 
tion from noise or they may be a population of transient objects in 
the ambient cloud, not related to the active star formation. 

The plots of the velocity dispersion against radius in Fig.[2T]do 
not show an obvious correlation. This is reflected in the low values 
of the correlation coefficients (0.2 for each cloud). In these plots the 
populations of low-mass condensations mentioned above are also 
plotted in red, but here they are not distinct from the other conden- 
sations. Several other surveys have also found very weak or no cor- 
relation between <7 V and R for condensations in molecular clouds 
(e.g. |Kramer et al.|19961|Onishi et al |2002] ). For a larg er ranse of 
molecular cloud sizes, |Larson| ( |1981| ) found that a v oc /r° - 4 . The 
lack of correlation found here could be because of the small range 
in R and ov and the large scatter in the values. 

Fig.[22]shows the relationship between M v i r and Mlte- Again, 
the low-mass condensations identified above are plotted in red, and 
they appear as a distinct population in this plot. The rest of the 
condensations are relatively well-correlated, with correlation coef- 
ficients of 0.9 for NGC 2024 and 0.8 for NGC 2071. The plots 
appear to be well-fitted with a power law (solid line), M v i r oc M^ TE , 
with p = 0.4 for NGC 2024 and 0.6 for NGC 2071, calculated 



& Sunada 



2009) found M vir /M LT E 1 



M- T E 33 for H 1J CO+ con- 
densations in Orion B (i.e. /3 = 0.67), which is very close to our 
value for NGC 2071. The dotted line on each plot shows where the 
condensations are in approximate equipartition {M N - n = Mlte)- The 
condensation virial masses are in fact all much higher than their 
LTE masses, implying that the condensations are unbound. How- 
ever, there are large uncertainties in both the virial and LTE masses 
(because of, for example, uncertainties in the distances to the cloud, 
the fractional abundance of 13 CO, the excitation temperature, and 
the assumed density profile), therefore we cannot make any definite 
conclusions about the fate of any particular condensation from this 
plot. However, the relative positions of the condensations to each 
other should be more robust and suggests that condensations from 
the low-mass population (in red) are in general less bound than the 
rest (i.e. they have higher M v i r /MLTE ratios). 



5.3 Condensation mass function 

The differential condensation (or core) mass function (CMF) is 
usually fitted by the power law 

dN 



dM 



oc M 



(20) 



If these condensations are the direct precursors of stars, the form 
of the CMF may provide insights into the origin of the stellar ini- 
tial mass function (IMF), which is described by Equation |20| with 
a = 2.35 ( Salpeter 1955 ) over a large range of environments. For 
instance the similarity of the CMF to the IMF for the earliest stages 
of protostellar evolution, may rule out models of star formation 
where the shape of the CMF is set at later stages (e.g. Bate & Bon- 
|nell|2005|L 



Fig. 23 compares the plots of dN /dM versus Mlte for both 
clouds. The crosses mark the mass function including all of the con- 
densations and the diamonds mark the mass function with the pop- 
ulation of low-mass condensations, discussed in Sec. |5.2| removed. 
The errors plotted show the statistical uncertainties of y/N, where 
N is the sample number in each mass bin. The two mass functions 
for each cloud are the same for Mlte ^ 0.02 M . Below this mass, 
the total mass function continues to increase with decreasing mass, 
whereas the mass function with the low-mass population removed 
has a turning point and begins to decrease with decreasing mass. 
Due to the increasing incompleteness in the mass function at these 
low masses arising from our detection thresholds, we have only at- 
tempted to fit the mass function for Mlte ^ 0.02 M . 

Reduced- squared fitting of a single power law mass function 
to NGC 2024 gives a ~ 1.3 for M LT e > 0.02 M (plotted with the 
dashed line on Figure [23}. The data suggest that there is a break 
in the power law at Mlte ~ 2.0 M . Fitting a double power law 
gives «i ow ~ 1.0 for 0.02 < M L te/M < 2.0 and «high ~ 2.6 for 
Mlte >2M . 

For NGC 2071, a single power law gives a ~ 1.7 for Mlte ^ 
0.02 M , but it is clear that the mass function flattens off for 0.02 < 
Mlte/M < 0.12. Fitting a double power law gives (X\ ow ~ 0.06 
for 0.02 < M L te/M < 0.12 and a high ~ 2.3 for M L te > 2M . 
The data suggest there could be another break in the power law at 
M L te ~ 1.0 M . 

The single power law values for a are similar to the values 
of 1.6 — 1.8 derived using isotopologues of CO by |Kramer et al.| 
(1998J for seven molecular clouds, including Orion B South. For 
the double power law fits, the values of «high are similar to those 
derived by J01 for the dust condensation of Orion B North and 
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Figure 21. Mlte — ^ (top panels) and o v —R (bottom panels) relations for 
NGC 2024 (left) and NGC 2071 (right). Condensations with M LTE /M < 
R/pc are plotted in red because they appear to form a separate population. 
In the top panels the dotted lines show M LTE /M = R/pc, the dashed lines 
show Mlte 06 R 2 (Larson's Law), and the solid lines show the lines of best 
fit: M LTE « R 2 6 for NGC 2024 and ~ R 1 7 for NGC 2071. 



Ikeda et al. 



|2009| > for the H 13 CO + condensations of Orion B, 



by 

which had values of ranging from 2.3 to 3.0. These values are also 
very similar to the slope of the initial mass function (a = 2.35), 
which has been interpreted as evidence for a physical link between 
the CMF and IMF . 

Iked a et aT] ([2009 ) also found (X\ ow ~ 0.06, just as we have 
found for NGC 2071 . They suggested the flattening could be a con- 
fusion effect where low-mass cores are misidentified as parts of 
high-mass cores. The turnover point found by Ik eda et ah] {2009) 
was much higher however, at 6 M , compared to our value of 
0.12 M©. Again this could be because the 13 CO LTE masses are 
underestimated due to opacity and freeze-out. 

From our analysis on the 13 CO condensations in NGC 2024 
and NGC 2071, we have found that NGC 2024 has a wider mass 
range than NGC 2071 (with masses up to ~ 13 M for NGC 
2024, compared to ~ 5 M for NGC 2071). Compared to the virial 
masses and condensation masses derived from other tracers, the 
13 CO mass estimates appear to be underestimated, possibly due 
to optical depth and freeze-out effects. The relationships between 
the condensation properties are similar to those seen in previous 
surveys of condensations of CO isotopologues andH 13 CO+. The 
slopes of the CMF at the high-mass ends (~ 2.6 for NGC 2024 and 
~ 2.3 for NGC 2071) are very similar to the slope of the initial 
mass function. 



6 SUMMARY 

We have offered a first look at the data on Orion B being collected 
by HARP/ACSIS on the JCMT for the Gould Belt Legacy Sur- 
vey. Our observations of CO, 13 CO and C 18 / = 3 -> 2 provide a 
comprehensive determination of the characteristic physical proper- 
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Figure 22. Plots of virial mass versus LTE mass for NGC 2024 (left) and 
NGC 2071 (right). Condensations with M L te/M < R/pc are plotted in red. 
The solid lines on each plot show the lines of best fit (M v j r oc M^ E for NGC 
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2024 and oc M^ E for NGC 2071) and the dotted lines show M y 
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Figure 23. The differential mass functions of 13 CO condensations in NGC 
2024 (left) and NGC 2071 (right). The crosses mark the mass function in- 
cluding all of the condensations and the diamonds mark the mass function 
with condensations with M L te/M < R/pc removed. The errors plotted 
show the statistical uncertainties of y/N, where N is the sample number in 
each mass bin. The dashed lines show the single power law fits (a — 1 .3 for 
NGC 2024 and a = 1 .7 for NGC 2071), and the solid lines show the double 
power law fits (cc low = 1.0 for 0.02 < M L te/M < 2.0 and a high = 2.6 for 
Mlte > 2M for NGC 2024 and ai ow = 0.06 for 0.02 < M LT e/M < 0.12 
and a high = 2.3 for M LT e > 2M for NGC 2071). 



ties and dynamics of star forming regions. We have observed two 
large regions in Orion B, (10.8x22.5) arcmin 2 in NGC 2024, and 
(13.5 x21.6) arcmin 2 in NGC 2071. Towards NGC 2024, 12 CO and 
13 CO are detected throughout the region, while C 18 is detected in 
an extended ridge, which follows the dust emission, and also in 
a fragmented ring surrounding a bright optical nebula. The 12 CO 
line profiles indicate multiple line of sight components. Compar- 
isons of the intensity ratios of the three isotopologues indicate that 
the 12 CO is optically thick throughout the cloud, and also in the 
lower- velocity outflow material, while the C 18 is generally opti- 
cally thin. These opacity results also apply to NGC 2071. 

Towards NGC 2071, the molecular gas is concentrated in a 
region almost completely surrounding a bright optical nebula. The 
12 CO emission is heavily self-absorbed, and 13 CO also shows self- 
absorption in the densest regions. The main outflow is very ener- 
getic, and is seen to produce extended wings even in the C 18 line 
profiles. 

Both clouds show a complex clumpy and filamentary struc- 
ture. In the regions traced by C 18 0, the equivalent widths suggest 
the bulk of the quiescent material has a lower velocity dispersion 
than the material in which it is embedded, increasing only in the 
regions of current star formation. The energetic outflow activity 
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is contributing a few percent of the kinetic energy of NGC 2024, 
while the outflows dominate the kinetic energy towards NGC 207 1 . 
Towards NGC 2024, the low-velocity red-shifted material is ex- 
tended and filamentary, while the low- velocity blue- shifted mate- 
rial is clumpy and compact. 

A CLUMPFIND analysis on the 13 CO data finds multiple 
13 CO condensations spatially associated with SCUBA dust cores, 
-7.4 towards NGC 2024, and -3.7 towards NGC 2071. The con- 
densations in NGC 2024 have higher LTE and virial masses than 
in NGC 2071, and towards both clouds, M v i r ^ Mlte- The slopes 
of the condensations' mass function at the high-mass ends are in 
agreement with the slope of the initial mass function. Towards NGC 
2024, we also detect a group of condensations that appear to be less 
bound, which may be tracing transient clumps that do not go on to 
form stars. 

These data are still being collected, and the Gould Belt Sur- 
vey has a final la noise requirement of of 0.3 K in 1.0 kms -1 for 
12 CO, and 0.3 K in 0.1 kms -1 for C 18 per 7.5 arcsec pixel. 13 CO 
is obtained simultaneously with C 18 0, and the C 18 requirement 
sets a la noise level of 0.25 K in 0.1 kms -1 for 13 CO. When the 
full HARP data set has been obtained, we will carry out more de- 
tailed analyses of the individual cores and outflow properties in the 
Orion B clouds. The completeness limit for the 13 CO results will be 
improved, although overall condensation masses, and the slope of 
the condensations' mass function will not be significantly affected. 
With datasets of the same sensitivity, the properties of the Orion 
B clouds will be compared to the other star-forming clouds being 
observed by the GBS. This data set will be highly complementary 
to the SCUBA-2 and POL-2 observations planned as part of this 
survey, and to planned Herschel key programmes. 
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